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Abstract 

We consider nonlinear, scaling-invariant = 1 boson+fermion supersymmetric 
systems whose right-hand sides are homogeneous differential polynomials and satisfy 
some natural assumptions. We select the super-systems that admit infinitely many 
higher symmetries generated by recursion operators; we further restrict ourselves 
to the case when the dilaton dimensions of the bosonic and fermionic super-fields 
coincide and the weight of the time is half the weight of the spatial variable. We 
^ . discover five systems that satisfy these assumptions; one system is transformed to 

the purely bosonic Burgers equation. We construct local, nilpotent, triangular, 
' weakly non-local, and super-recursion operators for their symmetry algebras. 
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I Introduction. We consider the problem of a complete description of = 1 nonlinear, 

scaling invariant evolutionary super-equations {ft = (p'^ , W = 0''} that admit infinitely 
many symmetries {fs = F, bg = B} proliferated by recursion operators TZ; here b is the 
bosonic super-field and / is the fermionic super-field. The axioms for selecting = 1 
nonlinear homogeneous polynomial evolutionary systems with higher symmetries were 
_c3 I suggested |^ by V. V. Sokolov and A. S. Sorin; the axioms are discussed in 

By construction, the equations are scaling invariant: their right-hand sides are differ- 
ential polynomials homogeneous w.r.t. a set of (half-)integer weights [6] = — |, [x] = —1, 
[t] < 0, [/], [b] > 0; we also assume that the negative weight [s] is (half-)integer. Here 
we denote by 9 the super- variable and we put Oq = Dg + 9 such that = D^, here 
Dg and are the total derivatives w.r.t. 9 and x, respectively. All notions and notation 
follow 0, see also |2I for details. 

In this paper, we investigate the properties of systems such that the weight of the time 
t is [t] = —J. We also assume [/] = [&] = | (the weights may not be uniquely defined). 

The first version of SsToOLS package for Reduce was used for finding the systems 
that satisfy the above axioms and possess higher symmetries under the bound —5 ^ 
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[s] ^ — |. Five systems were thus discovered, see Table I below. Later, we used the 
second version of SsToOLS [Hj for symmetry analysis of the super-systems in P| and for 
constructing conservation laws and recursion operators for their symmetry algebras. The 
method of Cartan forms j3] for the recursion operators was applied. Within this approach, 
the recursions are regarded as symmetries of the linearized equations. Namely, we 'forget' 
the internal structure of the symmetry flows fs = F, bg = B and operate with F and B 
as we do with the components of solutions of the linearized equations. The expressions 
TZ = R{F, B) are the recursion operators if each TZ satisfies the linearized equation again 
and if they are linear w.r.t. F, B, and their derivatives. 

Let us introduce some notation. Assume 7?. is a recursion for an equation and consider 
the symbol orcT'^^'^ieight • '^^^ subscripts 'ord' and 'weight' denote the differential order 
and the weight of the recursion 71, respectively, and the superscript 'layers' (if non-empty) 
indicates the required number of layers of the nonlocal variables assigned to conservation 
laws. The symbol 'jj' denotes the number of recursions for a given differential order, 
weight, and the nonlocalities. Further, we denote by L the local recursion operators, by 

the nonlocal or weakly non-local ^1 IH] recursions, the symbol Z denotes a nilpotent 
recursion whose powers equal zero except for a finite set, and is a super-recursion that 
swaps the parities of the flows. 

Now we list the flve new super-equations and indicate their recursions. The weights of 
the recursion operators are calculated w.r.t. the standard values [/] = [b] = |, [t] = — |. 
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bt = b^ + def 
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\ ft = deb + fb, 
bt = def 
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Table L 

It turns out that these equations exhibit practically the whole variety of properties 
that superPDE of mathematical physics possess. Let us discuss the properties of the 
equations present in Table I in more detail. 

1. The Burgers equation. First we construct an = 1 supersymmetric represen- 
tation of the Burgers equation and investigate its properties. We consider the system 

ft = deb, bt = b^ + def. (1) 

There is a unique set of weights [/] = [b] = ^, [t] — [x] = — 1 in this case. Hence 
we conclude that the above system precedes the invariance w.r.t. the translation along x. 
Equation (0) admits the continuous sequence Q of higher symmetries fs = (t)^,bs = 4>^ 
at all (half-) integer weights [s] < — |. Also, there is the continuous sequence (jU of 
supersymmetries for Eq. at all (half-) integer weights [s] < — | of the fermionic 'time' s. 

System (jT)) is obviously reduced to the purely bosonic Burgers equation b^ = btt — 2bbt. 
We emphasize that the role of the independent coordinates x and t is reversed w.r.t. the 
standard interpretation of t as the time and x as the spatial variable. The Cole-Hopf 



substitution b = —u^^ut from the heat equation Ux = Uu is thus the solution for the 
bosonic component of (^. 

Further, we introduce the bosonic nonlocahty w of weight [w\ = by the rules dew = 
— /, Wt = —b. The variable w is a. potential for both fields / and b. The nonlocality satisfies 
the potential Burgers equation = Wu + wf such that the formula w = hiu gives the 
solution; the relation / = —dgw determines the fermionic component in system (|T|). 

Now we extend the set of dependent variables /, b, and w by the symmetry generators 
F, B, and W that satisfy the linearized relations upon the flows of the initial super-fields, 
respectively. In this setting, we obtain the recursion 

^ _ fFx~dgfF + fxW\ fDx-dgf + fxd^' \ 

'^[^1 - -defB + bx W) ^ bx d^' Dx - def) 

of weight [sr] = — 1. In agreement with pp, the above recursion is weakly non-local We 
recall that a recursion operator R is weakly non-local if each nonlocality dg^ is preceded 
with a (shadow |3j of a nonlocal) symmetry and is followed by the gradient ipa of a 
conservation law: R = local part -|- J2a ■ o ijj^. From P it follows that this property 
is automatically satisfied by all recursion operators which are constructed using one layer 
of the nonlocal variables assigned to conservation laws. 

Recursion Q generates two sequences of higher symmetries for system (^: 



/A / debx - def deb - fxb \ f fx\ f fxx - '^deffx 

bj ^ \def. - {defy - b^def + bbx) ^ ' ' ' ' [bxj ^ [bxx - 2defbx 

Also, recursion (j2)) produces two infinite sequences of supersymmetries for 

'def\ (defx-{def?~U\ ( fdeb-bOef + bx 

deb) ^ [debx - def deb - bxf ) ^ ' ' ' ' \bdeb - f def + fx - fb^ 



(3) 



(4) 



Remark 1. System is not a supersymmetric extension of the Burgers equation; it 
is a supersymmetric representation of the Burgers equation. However, symmetries Q 
and (jH) are not reduced to the purely bosonic (x, t)-independent symmetries |3] of the 
Burgers equation (particularly, owing to the interchanged role of the variables x and t). 
We finally recall that the Burgers equation has infinitely many higher symmetries that 
depend explicitly on the base coordinates x, t but exceed the set of axioms [2; we use. 

Two supersymmetric generalizations (A^ = and = 2) of the Burgers equation are 
constructed in j2]. The A^ = extension relates it with integrable flows on associative 
algebras. The N = 2 Burgers equation contains a KdV-type component and admits an 
N = 2 modified KdV equation symmetry flow. 



2. A system with nonlocal recursions. The second system, 

ft = deb + fb, k = def, (5) 

is also homogeneous w.r.t. a unique set of weights [/] = [b] = |, [t] = [x] = — 1. 
Similarly to the supersymmetric representation for the Burgers equation, Eq. (jSj) 
admits symmetries {fs, bs) for all weights [s] < — |. 

We conjecture that system (0) has only one conservation law that defines the fermionic 
variable w of weight hj Wt = f, dew = b. Then, many nonlocal conservation laws and 



hence many new variables appear. We use the fermionic variable v whose weight [v] = | 
is minimal: we set Vt = deb ■ wfb + fxwf and dgv = —dgb ■ fb + dgf ■ dgb ■ w + b^wf. Now, 
there are nontrivial solutions to the determining equations for recursion operators. First, 
we obtain the recursion of zero differential order with nonlocal coefficients: 

_ /-deb ■ wfB + wvF + v ■ B\ 
[-i^l ~ V debwfF -vF + VW B )' 

Also, we get a nonlocal operator with nonlocal coefficients, 

^ (deb Vw - def deB wf - def deb Wf + def dgb Fw + defV+ Vwfb\ 
f"^' V deBdebwf + debV-debFwfb + defdebVwf + defVw-Vfb J' 

The coefficients of the recursions found for [s/j] = — 2| and [sr] = —3 are also nonlocal. 



3. A triplet of super-systems. Finally, we consider the three systems 

ft = -afb, k = b^ + def (6) 

which differ by the values a = 1, 2, and 4 of the parameter a and therefore exhibit rather 
different properties. The weights for the above equation are multiply defined, and we 
choose the tuple [/] = [b] = |, [t] = — |, [x] = — 1 to be the primary 'reference system.' 

Case a = 2. First, we fix a = 2 and consider Eq. (0): we get ft = —2fb, bf = 
b^ + def. The weights for symmetries are [s] = — |, [s] = —1, and then Eq. (0) admits a 
continuous chain of symmetry flows for all (half-)integer weights [s] < — 2|. Surprisingly, 
no nonlocalities are needed to construct the recursion operators, although there are many 
conservation laws for this system. We obtain purely local recursion operators TZ that 
proliferate the symmetries: <f = {F, B) \^ Lp' = TZ for any (p. The recursion 

/ ^a,F def f + lldeF fb^ + l^deffF + ?>def Fb'' + \fxFf \ 

''^[-^1 " lldeB fb' + 8deb Fb' + 22deb fBb + 7{deffB + 

\ Udef Bb' + fdef deB f + \def deb F + \bxF f + f^Fb + hfJB j 

of weight [sb\ = —2 is triangular since W does not contain B. Also, we obtain the 
recursion of weight 2^; its components are 

1 = -2debFfb' - deFdef fb - deF fb^ - \f,Ffb - 2def fBb', 

7^['_^lJ = deBfb' + debFb' + deb fBb' + ^def.Ff + 

' + IdeFb' + \deF{def)' + deFdefb' + ^deF fj + {deffBb + 
+ def Bb^ + def deB fb + def deb Fb + def deb fB + pef FJ + 
+ Idef f.F + \bxFfb + \FJb' + \fxFb' + \fJBb. 

Further, we get a triangular nilpotent operator of weight —3 such that "^[^3] = and 
^[_3] = {defyfF + 6{def)'fb'F+12deffb'^F + 8fb^F. The above recursion is a re- 
currence relation (21 which is well-defined for all symmetries of Eq. ©. Another local 
recursion for [s] = —3 is huge and therefore omitted. 



For a = 2, system ^ admits at least three super-recursions ^{R^ , B!') such that the 
parities of and R^ are opposite to the odd parity for / (and hence for F) and to 
the even parity of b and B. This property is possible owing to the presence of the odd 
variable sr. The triangular zero-order super-recursions are T^f_2] = ^dofFfb + SFfb^, 
= -AdebFfb + 2{deffF + 6def Fb"^ + Adef fBb - f^Ff + AFb'^ + 8fBb^ and 

^ ^( ~doff.F-2f,Fb' \ 

\deb f,F - def b,F + def f,B - 2b,Fb^ + 2f,Bby 

for weights [s/jj = —2 and [s/jj = — 2^, respectively; the third super- recursion found for 
[s/j] = — 2| is very large. Quite naturally, system (jH} has infinitely many supersymmetries 
if a = 2. 



Case a = 1. For a = 1 from (jHl) we obtain the system ft = —fb, bt = b"^ + def . The 
default set of weights is the same as above: [/] = [^] = |, [i] = — |, and [x] = —1. The 
sequence of symmetries is not continuous and starts later than for the chain in the case 
a = 2. We find out that there are symmetry flows if either [s] = [t] = (the equation 
itself), [s] = [x] = —1 (the translation along x), or [s] < — 3| such that a continuous chain 
starts for all (half-)integer weights [s]. 

Similarly to the previous case, no nonlocahties are needed to construct the recur- 
sions, which therefore are purely local. The recursion operator TZ-^ = 0, TZ^ ^ = 

(deff Ff + 3dgf Ffb"^ + ^Ffb^ of maximal weight [sr] = -2\ is nilpotent: = 0. For 
the succeeding weight [sr] = —3, we obtain a nilpotent local recursion with components 

7^f_3, = IdeFidofff + IdgFdgffb' - lidgffF - l{dgf)'Fb' + 

+ hdef deb Ffb + fa,/ f^Ff + f /.F/fe^, 
7^|'_3] = def.Ffb-^deFdebfb^ - lfdeFdefdebf + ndeFfJb + 
+ lideffdeBf + l{dgf)'debF + Idef deB fb' + IdefdebFb' - 
- 559,/ deb fBb + fdef b.,F f + def fJB + fb^Ffb' + ^fJBb\ 

It generates symmetries of system ©; the differential order of 7?.[_3] is positive. 

Case a = 4. Finally, let a = 4; then system © acquires the form ft = —4/6, 
bt = b'^ + def. Again, the basic set of weights is [/] = [b] = i, [t] = — i, [x] = —1, and 
system © admits the symmetries {fs, bs) such that their weights are [s] = —1 or 
[s] < — 3| w.r.t. the basic set. This situation coincides with the case a = 1. Again, no 
nonlocahties are needed for constructing the recursion of minimal weight [sr] = — 3|: 

, = -UdebFfb' - deF {defffb - AdeFdeffb' - 3deF fb' - 
- Aide f)' fBb' - AdefdebFfb' - peff.Ffb-UdeffBb' - 2f,Ffb^ 
7^|' = SdeBfb"^ + SdebFb' + ddgbfBb' + \def. defFf - \def. Ffb' + fa.Ffoe + 
' + deF deb fb' + \deF {deff + \deF {deffb' + IdgFdgfb^ + deF def deb fb + 
+ ^^dgF def fj + \deF fjb' + (def fBb + A{dgf)'Bb' + {dgffdgB fb + 

+ [dgffdebFb + {deffdgbfB + l{dgffFJ + \{dgfff.,F + ?,dgfBb'> + 

+ Adef deB fb"" + Adef deb Fb' + lOdef deb fBb' + ^def b^Ffb + dgf FJb' + 
+ Idgf f,Fb' + Idgf fJBb + 2b^Ffb' + FJb^ + \f,Fb^ + fJBb\ 



No nilpotent recursion operators were found for system if a = 4. 

Remark 2 . We discovered that an essential part of recursion operators for supersymmetric 
PDE are nilpotent. At present, it is not clear how the nilpotent recursion operators con- 
tribute to the integrability of supersymmetric systems and what invariants they describe 
or symptomize. Further, we emphasize that this property does not always originate from 
the rule '/ ■ / = 0', but this is an immanent feature of the symmetry algebras. More 
generally, the nilpotent recursions are quite natural in the bosonic sector, too. We have 

Example (I. S. Krasil'shchik, private communication). Consider a system of linear or- 
dinary differential equations x = A{t) x. Then any nilpotent constant matrix R that 
commutes with the matrix A is a recursion. 

It would be of interest to construct an equation S that admits nilpotent differential 
recursion operators {Ri, . . . \ R^' =0} which generate an infinite sequence of symmetries 
ip, Ri^ o Ri^(ip), ... for S. Here we assume that at least two operators (without 

loss of generality, Ri and R2) do not commute and hence the flows never become zero. 
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